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We study the redshift drift, i.e., the time derivative of the cosmological redshift 
in the Lemaitre-Tolman-Bondi (LTB) solution in which the observer is assumed to 
be located at the symmetry center. This solution has often been studied as an anti- 
Copernican universe model to explain the acceleration of cosmic volume expansion 
without introducing the concept of dark energy. One of decisive differences between 
LTB universe models and Copernican universe models with dark energy is believed 
to be the redshift drift. The redshift drift is negative in all known LTB universe 
models, whereas it is positive in the redshift domain z < 2 in Copernican models 
with dark energy. However, there have been no detailed studies on this subject. In 
the present paper, we prove that the redshift drift of an off-center source is always 
negative in the case of LTB void models. We also show that the redshift drift can 
be positive with an extremely large hump-type inhomogeneity. Our results suggest 
that we can determine whether we live near the center of a large void without dark 
energy by observing the redshift drift. 
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I. INTRODUCTION 

The standard cosmological model is based on the so-called Copernican principle that 
we are not located in a special position in the universe. This model can naturally explain 
almost all observational data, and consequently seems to imply that the Copernican principle 
is a reality. However, we should not blindly rely on this principle without observational 
justifications. 

The highly isotropic cosmic microwave background (CMB) implies the isotropy of our 
universe at our position. Then, if we assume that the Copernican principle is correct, 
we necessarily arrive at the conclusion that our universe is isotropic at every position, or 
equivalently, our universe is homogeneous and isotropic. By virtue of the homogeneity and 
isotropy, the standard cosmological model is determined by several parameters called the 
cosmological parameters. In order to determine the cosmological parameters, observational 
data are interpreted under the assumption of a homogeneous and isotropic universe on 
average. Here, we should note that it is not clear at all how large the systematic errors 
would be in the determination of the cosmological parameters, if the Copernican principle 
is abandoned. Thus, it is an unavoidable task in observational cosmology to investigate 
possible "anti-Copernican" universe models and test if such models can be observationally 
excluded. 

Almost all anti-Copernican universe models are based on the Lemaitre-Tolman-Bondi 
(LTB) solution which describes the dynamics of a spherically symmetric dust. In anti- 
Copernican universe models, an observer like us is usually assumed to be located in the 
neighborhood of the symmetry center [l -^|. In recent years, such models have attracted 
much attention, since some LTB universe models can recover the observed distance-redshift 



relation without introducing dark energy |6H22j, and various ways to observationally test 



23|-|58|. 



these models have been proposed by many authors 

As shown by previous studies, the LTB solution can explain various observational data 
besides the distance-redshift relation without the need to introduce dark energy. This is 
because this solution has functional degrees of freedom with respect to the comoving radial 
coordinate. In order to check the LTB universe models observationally, it is crucial to find 
observable quantities which can reveal differences between the LTB universe models and 
Copernican universe models with the dark energy. One such quantity is believed to be the 
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redshift drift, i.e., the time derivative of the cosmological redshift 46|. In the case of the 
ACDM model, which is the most likely Copernican model at present, the redshift drift is 
positive in the redshift domain z < 2, since the cosmological constant A causes repulsive 
gravity. By contrast, there is no exotic matter with the violation of the strong energy 
condition in the LTB solution. Thus, as long as there is no highly inhomogeneous structure, 
the redshift drift might be negative in LTB universe models. Although several authors have 



pointed out the importance of the redshift drift |41|, |46|, |59|, |60| , there has been no detailed 
study of its general behavior in LTB universe models. It is the purpose of this paper to 
investigate it. 

This paper is organized as follows. In Sec. HIl we briefly review the LTB solution. In 
Sec. IHIl we derive the equation for the redshift drift and show the behaviour of the redshift 
drift near an observer located at the symmetry center in LTB universe models. In Sec. IIV| 
we define LTB void models and prove a theorem on the redshift drift in these models. In 
Sec. |Vl we show that the redshift drift can be positive even in an LTB universe model, if 
an extremely large hump-type mass density distribution exists. Sec. |VT] is devoted to the 
summary and discussion. 

In this paper, we denote the speed of light and Newton's gravitational constant by c and 
G, respectively. 



II. THE LTB SOLUTION 



As mentioned in the introduction, we consider a spherically symmetric inhomogeneous 
universe filled with dust. This universe is described by an exact solution of the Einstein 
equations, which is known as the Lemaitre-Tolman-Bondi (LTB) solution. The metric of 
the LTB solution is given by 

ds > = _ c 2 rft 2 + (d r R(t,r)f ^ + r2 2 
1 — k[r)r l 

where k(r) is an arbitrary function of the radial coordinate r. The matter is dust whose 
stress-energy tensor is given by 

TV = pu^u", (2) 

where p = p{t,r) is the mass density, and w M is the four-velocity of the fluid element. The 
coordinate system in Eq. ([!]) is chosen in such a way that = (c, 0, 0, 0). 
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The circumferential radius R(t, r) is determined by one of the Einstein equations, 

where M(r) is an arbitrary function related to the mass density p by 

/ x I dM 

P ^ r) = 4^R^- (4) 

M(r) is known as the Misner-Sharp mass that is the quasi-local mass naturally introduced 
into the spherically symmetric spacetime {(51 1 . Here, it should be noted that the Misner- 
Sharp mass is not necessarily a non- decreasing function of the comoving radial coordinate 
r, even if the mass density p is non-negative. In the case that a hypersurface labeled by t 
is a homogeneous and isotropic space with positive curvature, r can be chosen so that the 
circumferential radius is given by R = a(t) sinr, where a(t) is a positive function of time (in 
this case, k(r) = r~ 2 sin 2 r). Thus, d r R = a(t)cosr is positive for < r < tt/2, whereas it 
is negative tc/2 < r < n. Since p is spatially constant by assumption, the derivative of the 
Misner-Sharp mass dM/dr is positive for < r < tt/2, whereas it is negative for tt/2 < r < tt. 
However, in this paper, we assume that the Misner-Sharp mass is a monotonically increasing 
function of r in the domain of interest. This assumption is equivalent to the one that d r R 
is positive if p is positive. 



Following Ref. 



62| . we write the solution of Eq. ([3]) in the form, 



R(t,r) = (6GMy> s \t-t B (r)] 2 > 3 S(x), (5) 

where £b(?") is an arbitrary function which determines the big bang time, and S(x) is a 
function defined implicitly as 



S(x) 



cosh J— 77 — 1 — (sinh J—rt — \/—ri) 2 ^ 3 
X ! ■ x = : X ! X f or x < o 

6 1 / 3 (sinh ^—f] — \/—i]) 2 / 3 ' 6 2 / 3 
1 - cos y/rj _ {^n-sm^qfl 3 

6 1 /3( v ^_ siriv ^)2/3 ' X 6 2/3 IOT X>U - 

(7) 

The function S(x) is analytic for x < (7r/3) 2//3 . Some characteristics of the function S(x) 
are given in Refs. 60|] and 62]. 

As shown in the above, the LTB solution has three arbitrary functions, k(r), M(r) and 
£b(?")- One of them is a gauge degree of freedom for the rescaling of r. In this paper, since 
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M is assumed to be a monotonically increasing function of r, we can fix this freedom by 
setting 

M = ^ W 3 , (8) 

where po is the mass density at the symmetry center at the present time to, i.e., po = p(to, 0). 
By this choice, we have 

R(t ,r)=r + O(r 2 ). (9) 

As in the case of the homogeneous and isotropic universe, the present Hubble parameter H 
is related to p as 

Hi + k(0)c 2 = ^vrGpo. (10) 

III. EQUATION FOR THE REDSHIFT DRIFT 

In order to study the cosmological redshift and the redshift drift, we consider ingoing 
radial null geodesies. The cosmological redshift z of a light ray from a comoving source at 
r to the observer at the symmetry center r = is defined by 

{ >- *'(A(0)) ' 1 ' 

where k l is the time component of the null geodesic tangent, and A is the affine parameter 
which can be regarded as a function of r. From the geodesic equations, we have the equation 
for the redshift z as 

ch = (1 + z)d t d r R 
dr ~ cVl - kr 2 ' 

The null condition leads to 

dt d r R , . 

-r = ; !3 

dr cy/l - kr 2 V ' 

We denote the trajectories of light rays observed by the central observer at t — t and 

t = t + 5t , respectively, by 

z = z lc (r;t ) ^ 
^ t = ti c (r;t ) 
and 

z = zi c (r; t + 5t ) =: zi c (r; t ) + Sz(r) 

(15) 

t = ti c (r; t + St ) =: t k (r; t ) + St(r) 
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Here, by their definitions, we have t\ c (0;t ) = t , zi c (0;t) = 0, 8z(0) = and 5t(0) = St . 
Substituting Eq. f lT5l) into Eqs. f|T2|) and ( JTBl . and regarding 5z{r) and St(r) as infinitesimal 
quantities, we obtain 

i = BAR (1+06^ 
dr C y/l - kr 2 cy/l - kr 2 

4-H = -=^£=St, (17) 
dr C y/i - kr 2 V ' 

where we have used the fact that ( THl) satisfies Eqs. ( |T2|) and ( fT3l) . and the arguments of 
d t d r R and d 2 d r R are i = ti c (r, t ) and r. 

Hereafter, we consider the case where the cosmological redshift z is monotonically in- 
creasing with r. We say that such a model is z-normal. Then, we replace the independent 
variable r by z = zi c (r;t ). By using 

d dz d (1 + z)d r d r R d 
dr dr dz cy/l — kr 2 dz ' 

we have 



(18) 



d . 5z d?d r R . 
dz 1 + z o t o r R 

— <yt = . 20 

rfz 1 + z v ; 



We can easily integrate Eq. ( 1201) to obtain 



ft = -^2-. (21) 
1 + z V ; 



By using the above result, Eq. (fl9l) is rewritten in the following form 

± / _ 1 d 2 d r R 

dz\l + z)~ {l + z) 2 d t d r R °" 1 j 

Here, let us study the redshift drift <5z in the neighborhood of the symmetry center. By 
the regularity at r = 0, we have 

k{r) = k + O(r), (23) 
t B (r) = <D(r), (24) 



where, by using the freedom of the constant time translation, we have set £b(0) = 0. From 
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Eq. d3J), we have 













t=t ,r=0 


d t d r R 


t=t ,r=0 


d 2 t d T R 


t=t ,r=0 


© and fi m0 = 



2 87rGpo „2 
-c fc H — = # ; 



47rGp 



O f7 2 



(25) 

(26) 
(27) 



from Eq.fl22J that 
d 



dz 



Sz 



d 2 t d r R 



t=to ,r=0 



dfd r R 



t=t ,r=0 



6t + O(z) = --n m0 6t + O(z). 



Thus, we have 



5z 1 

TT = --n m0 z + O(z 2 ). 

OlQ I 



(28) 



(29) 



The above equation shows that the redshift drift is non-positive near the center. This is the 
same behavior as that of the homogeneous and isotropic universe filled with dust. 



IV. THE REDSHIFT DRIFT IN LTB VOID MODELS 

We call an LTB universe model the LTB void model, if the following three conditions are 
satisfied. 

1. the mass density is non-negative; 

2. the mass density is increasing with r increasing in the domain r > on a spacelike 
hypersurface of constant t; 

3. d r R is positive; 

4. z- normality. 



Proposition 1 In LTB void models, d 2 d r R is negative. 



Proof. By Eq. ([3]), we obtain 



9 R ( l" r \ 

47lG JV V + 2 J P( f ^ X ) R2 (^ X ) d r R (^ X ) dx U ( 30 ) 
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where we have used Eq. (j4j) in the second equality. Since d r R is positive by the definition 
of LTB void models, we may replace the integration variable x by R = R(t, x) and obtain 

d R ( f R{t ' r) 
d 2 d r R(t, r) = 4vrG-^- - P R 3 + 2 P R 2 dR 

,d r R f RM fdp, , 



= - A * G ivl {M R+pR ) dR - (31) 

Since dp/dR = (d r R)~ 1 d r p is positive in the domain of R > 0, the integrand in the last 
equality of the above equation is positive. Q.E.D. 

Theorem In LTB void models, the redshift drift of an off-center source observed at the 
symmetry center is negative. 

Proof. Since the cosmological redshift z vanishes at r = 0, z is non-negative by the 
assumption of ^-normality. Further, the ^-normality leads to dtd r R > through Eq. (fl2|) . 
Then, since Sto > 0, we see from Eq. fl22|) that Proposition 1 leads to the following inequality 

d ( 5z 



, < 0- (32) 
dz \1 + z/ 

Since 5z should vanish at z — 0, we have 5z < for z > from the above inequality. Q.E.D. 



V. REDSHIFT DRIFT IN LTB UNIVERSE MODELS WITH A LARGE HUMP 

In the preceding section, we showed that the redshift drift observed at the symmetry 
center is negative for r > in LTB void models. Conversely, if there is a domain in which 
the mass density is decreasing with increasing r, the redshift drift might be negative. In this 
section, we show that it is true with hump-type mass density distributions. We consider the 
following two LTB universe models, 

(i) k(r) = and tB( r ) = /( r ; a , r i, r 2) 

with a = -1.7H \ n = O.UcHq 1 and r 2 = 0.9cH \ 

(ii) £b(?") = and k{r) = f(r;a,ri,r 2 ) 

with a = -100c- 2 # 2 , n = Olci^ -1 and r 2 = 0.2cH \ 
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where 



f{r;a,r 1 ,r 2 ) 



a(r — t*i) 3 (r\ — 57*ir 2 + lOr^ + 37*ir — 15r 2 r + 6r 2 
(r 2 - ri) 5 



for r < 7*1 , 
for 7*1 < r < r 2 
for 7*2 < r . 



(33) 

In FigsJU and EJ we show the redshift drifts of these models. Although we do not show the 
energy densities of these models, a large hump in the mass density distribution exists in 
each model as well as in £b(?") or k(r). Although there is a redshift domain with positive 
redshift drift in each example, the distance-redshift relations of these models do not agree 
with the observational data, and further, the inhomogeneities need to be very large. 





FIG. 1: The LTB universe model with k(r) = 0. The right panel depicts the redshift drift 5z/5to 
of the LTB universe model with k(r) = as a function of the redshift z. The other arbitrary 
function £b(?*) is shown in the left panel. 



VI. SUMMARY AND DISCUSSION 

In this paper, we studied the redshift drift in LTB universe models in which the observer 
is located at the symmetry center. We showed that, assuming that the mass density of the 
dust is positive, the redshift drift of an off-center source is negative if the mass density and 
the circumferential radius are increasing functions of the comoving radial coordinate. We 
also showed that if there is a very large hump structure around the symmetry center, the 
redshift drift can be positive. As a result, by observation of the redshift drift, we get a 
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FIG. 2: The right panel depicts the redshift drift 5z/5to of the LTB universe model with t^ir) = 
as a function of the redshift z. The other arbitrary function k(r) is shown in the left panel. 



strong constraint on void-type universe models: if the redshift drift turns out to be positive 
in some redshift domain, LTB void models can be rejected. 

By projects such as the Cosmic Dynamics Experiment (CODEX) 63M65I] in the Euro- 
pean Extremely Large Telescope (E-ELT) pa] and PREcision Super Stable Observations 
(ESPRESSO) in the Very Large Telescope array(VLT) [64], [65|, [67|, it is possible to get 
highly accurate spectroscopic observational data. The observability of the redshift drift by 



CODEX has been analyzed in Refs. 



64 



65. 



69(, and the test of LTB universe models 



6t 

by this project has been investigated in Ref. 41|. Targets of these projects are QSOs whose 



redshifts are larger than two. On the other hand, observing the sign of the redshift drift 
in the redshift domain z < 2 is very crucial to test LTB void universe models in contrast 
with Copernican universe models. This is because known examples of Copernican universe 
models with dark energy or modified gravity predict a positive redshift drift in the redshift 
domain z < 2, while it is negative for LTB void models. From, for example, the observation 



of compact binary stars by DECIGO 



7JH72J or BBO [73j, l74|, the redshift drift at z ~ 1 



can be measured 70J]. Observations by DECIGO or BBO over several years will make it 
possible to observe the sign of the redshift drift for z < 2 and to test LTB void models [75 ]. 
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